Hubble Constant, Lensing, and Time Delay in TeVeS 
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Hubble constant can be determined by the time delay of gravitational lensing systems. As data 
on time delay observation accumulates, it is time to revisit this approach. As in other dynamical 
phenomena in scales of galaxy and cluster of galaxies, gravitational lensing in these scales is also 
plagued by the problem of excess acceleration or gravity (a.k.a. missing mass problem). There are 
always some accelerations unaccounted for by luminous matter. Usually dark matter is introduced 
to interpret the discrepancy. However, Modified Newtonian Dynamics (MOND) is more successful 
in explaining the excess accelerations in galaxy scale. We adopt TeVe S as the relativistic version of 
MOND to study gravitational lensing phenomena, and we can evaluate the Hubble constant from 
the derived time-delay formula. To apply our method, we rely on the CASTLE quasar lensing 
survey and the Sloan Digital Sky Survey (SDSS). Four samples are suitable for our study. Using 
only the luminous part of the lensing galaxies, the average of the derived Hubble constant is 68.5 
km s _1 Mpc _1 . 



Hubble constant Hq is a long debated quantity in cos- 
mology for more than half a century. Basically, it comes 
from the relation between cosmological distance and re- 
ceding velocity of galaxies v = HqcI. Its inverse represent 
the age of the universe. 

The value of Hq is sensitive to the way we estimate 
the distance. Its value has been estimated by many 
distance-determination methods, such as Cepheids, tip 
of the red giant branch, maser galaxies, surface bright- 
ness fluctuations, Tully-Fisher relation, Type la super- 
novae, and other methods such as gravitational time de- 
lay, Sunyaev-Zel'dovich (SZ) effect, cosmic microwave 
background (CMB), see ref [HQ for details. In this work, 
we focus on gravitational lensing and time delay. Ob- 
taining Hq by Gravitational time delay was introduced 
in 1964 by Refsdal Q. Bright variable sources arc needed 
and Refsdal suggested supernovae. Subsequent works 
on Hq, however, used quasar lensing. Now there are 
19 systems of quasar lensing with time delay measure- 
ment dij. 

One advantage of using time delay to derive Hubble 
constant is that it is less sensitive to cosmological mod- 
els. Thus it provides a more direct probe to cosmological 
distance jlj. However, there are some uncertainties in 
determining mass distribution by image deflections and 
distortions from gravitational lensing. This is commonly 
known as "mass sheet degeneracy" [6J. Another source 
of uncertainty in mass is, of course, the missing mass 
problem. Missing mass is a long standing issue. Oort 
in 1932 and Zwicky in 1933 were the first to put forward 
the notion of missing mass in galactic systems @, § ■ The 
missing mass problem was neatly confirmed by the ob- 
served flat rotation curve in spiral galaxies [9h11| . Nowa- 
days, missing mass exists in nearly all galactic systems, 
clusters of galaxies, large scale structure and CMB. In 



fact, the problem should be interpreted in terms of ex- 
cess acceleration or gravity, i.e., there are some accel- 
eration cannot be accounted for by luminous matter. 
To compensate the excess acceleration one can, on the 
one hand, introduce dark matter into the system. On 
the other hand, one can modify Newton's law of mo- 
tion or law of gravity. In 1983, Milgrom proposed Mod- 
ified Newtonian Dynamics (MOND) to explain both flat 
rotation curve and Tully-Fisher relation [T3|. MOND 
asserts that when acceleration is smaller than about 
do = 10~ 10 m s~ 2 , Newton's second law of motion has to 
be changed: /2(|a|/oo)a = — V$n = »n, with fj,(x) w x 
for x -C 1, ft(x) — > 1 for x S> 1, where $n is the New- 
tonian gravitational potential. fl{x) is called the inter- 
polation function. For later calculation it is useful to 
introduce the inverted interpolation function v such that 
a = -V$ = i>(|a N |/a )a N 

MOND is very successful in explaining the dynamics 
of galactic systems [l3| . Recently, McGaugh showed that 
MOND can perfectly explain the Tully-Fisher relation in 
gas rich spiral galaxies without uncertainty parameter 
such as mass-to-light ratio [l4l [ToT ] . As usual the result 
created some debate Til- 1 8| . Nevertheless . many con- 
sider MOND is not quite successful on cluster of galaxies 
scale A recent study on the gravitational rcdshift 

of clusters of galaxies |20( has generated some debate 
on whether MOND is applicable to cluster scale, and it 
seems that MOND does not have difficulty in interpret- 
ing the data [HI]. In an y case, the original MOND is a 
non-relativistic theory and can not be applied to relativis- 
tic phenomena such as gravitational lens and cosmology. 
Two decades after the original paper by Milgrom [12j, 
in 2004 Bekenstein proposed the Tensor- Vector-Scalar 
(TeVe S) covariant relativistic gravity theory with MON- 
Dian dynamics as its non-relativistic limit |22j. In 2006, 
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Chiu et al. started from TeVe S and derived the cor- 
responding strong lens equation (23|. Later in 2009, 
Milgrom proposed another relativistic version of MOND 
called BiMOND El . It turns out that TeVe S and Bi- 



MOND have identical gravitational lensing results. The 
lensing theory has been applied to some galaxy lensing 
data, in which the galactic mass has been calculated and 
compared with population synthesis [25W27I . Once again, 
MOND can explain the data well. Moreover, the MON- 
Dian lensing mass is consistent with the dynamical mass 
from fundamental plane (28j . 

The observational tests above on galactic mass depend 
on the mass-to-light ratio, which cannot be measured di- 
rectly. Here we put forward another test, namely, Hubble 
constant Hq from time delay measurement. If MOND 
cannot get a reasonable value of Hq, it should be ruled 
out without any ambiguity. We solve the lens and time 
delay equation simultaneously to obtain both mass of the 
lens and the Hubble constant. Unlike dark matter the- 
ory with adjustable mass model, the mass distribution in 
MOND is fixed by the density profile of luminous matter. 

The discussion on light deflection due to gravitational 
force can be traced to Newton's Opticks. Modern view 
of light deflection is a relativistic gravitational effect, in 
which both time-like and space-like part in the metric 
contribute to the deflection angle. Using General theory 
of Relativity (GR) , Einstein derived the deflection angle 
which is just twice the Newtonian one. It turns out that 
the deflection angle derived from TeVe S or BiMOND 
is also twice the Newtonian one (i.e., the same as GR) 
except that the Newtonian potential is replaced by the 
MONDian potential. The deflection angle by a spherical 
lens in small angle approximation can be written as [23 . 
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where c is the speed of light, is the MONDian po- 
tential, go ps D\J) is the closest approach from the center 
of the lens, 9 is the image position, and g 2 = g 2 , + ( 2 . 
Z?l, D' L and Z?ls ar e the angular distances of the lens 
from the observer, observer from the lens and the source 
from the lens, respectively. The direction of the image 
is in the direction of the closest approach (projected on 
the sky). In thin lens approximation the corresponding 
lensing equation is given by 



P = 9+- a{9+) = a(6L) - 



(2) 



where /3 is the source position, 9± are the two images on 
the same and opposite side as the source, respectively, 
and a(9) = AipD-^s/Ds is commonly called the reduced 
deflection angle. 

Time delay is defined as the difference in time traveled 
by light along the actual path and along the undeflcctcd 
path. It can be derived from Format principle or from 



the geodesic equation in relativistic gravitation theory. 
As in deflection angle, the form of the time delay is the 
same for GR and MOND (with Newtonian potential for 
GR and MONDian potential for MOND), 
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The first and second term in Eq. ([3]) are referred to 
the geometric and the potential time delay. The angu- 
lar distance Da is proportional to the Hubble constant, 
Da = (c/Hq)Da{z), and Da depends on redshift z only. 
In the cases where the difference in time delay of two im- 
ages is available, the value of Hq (and the mass of the lens 
as well) can be obtained by solving difference time delay 
equation At = t(9^) — t(9 + ) and lens equation Eq. ([2]). 

It is interesting to note that the time delay difference 
of MOND for point mass lens in Bekenstein form ( = 
V<f>Ar(l + yj ao/| V$at| )) has exactly the same form as 
GR 



gpAt = 1 
D L D S (1 + z O 2 



^1 log 



(4) 



Here 9 E = ^AGMD LS /c 2 D h D s is called the Einstein 
radius. The only difference is that lens equation gives 
= 1 in GR, but in MOND it gives 



9 2 
7 E 



= 1, 



(5) 



where 9 = ^GM/aoD 2 (M is the mass of the lumi- 
nous matter). We can use 9o as a criterion for deep 
MOND regime if the gravitational acceleration at the 
closest approach is much smaller than ao, i.e., 9 3> 9q. 
Because in Eq. (0, 0| -> 0+0_ as 9± < 9 (Newto- 
nian) and 0| O (0+ + 0_)/tt as 0± > 9 (MONDian) 
, from Eq. (|4]) the value of H from MOND is smaller 
than that from GR, see Figure [T] In deep MOND regime 
the interpolation function is p,(a/ao) ~ a/do- If the ex- 
tend of the luminous matter is also much smaller than 
9 (i.e., practically a point mass), then the time delay 
difference is solely determined by the potential time de- 
lay, because the deflection angle approaches a constant 
in deep MOND regime [23j. The time delay difference 
is given by a simple formula which is independent of the 
choice of interpolation function, 



Dls H At = 1 



(G) 



We emphasize that Eq. ((6]) do not have any free param- 
eters for the interpolation function or mass of the lens, 
not even ao (the utmost important constant of MOND). 
ao remains in lensing Eq. ([5]) with mass of the lens only. 
Thus, if such a clean observation is available, then it has 
the potential to falsify MOND unambiguously. 
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In GR for Eq. ((4]), the first term is the geometric time 
delay and the second is the potential time delay. How- 
ever, in MOND the first term in Eq. (fj| comes from pri- 
mary MONDian potential time delay 7T0|(0 + - 0_)/20 o . 
Becuase this primary MONDian contribution cancel the 
geometric time delay, the remain is the first term in (fj| . 
When we go from Newtonian regime (8q/9± S> 1) to the 
deep MOND regime (0o/0± <C 1), 0e changes from the 
order of to the order of J 8q9±. The primary MON- 
Dian potential time delay remain finite, while the sec- 
ondary MONDian potential time delay (the second term 
in Eq. diminish. As a result, in the deep MOND 
regime the time delay (Eq. ©) becomes solely potential 
time delay. 



FIG. 1. Hubble constant (km s _1 Mpc -1 ) is plotted as a 

function of time delay difference (days) in a point mass lens model. 
The plot shows both the result of GR and deep MOND. Dashed 
line is GR, dotted-line is deep MOND in ACDM cosmology, solid 
line is deep MOND in massive neutrinos cosmology, and dot-dashed 
line for Hernquist model in MOND. To put the diagram into per- 
spective, we consider a quasar lensing object FBQ J0951+2635. Its 
measured time delay difference is 16 days and is marked by the ver- 
tical black line. For reference Hq = 72 km s _1 Mpc -1 is marked 
by the horizontal black line. 

The gravitational lens with time delay measurement 
that we are able to collect are not in the deep MOND 
regime. The gravitational acceleration ranges from 10~ 9 
ms -2 to 10~ 10 ms" 2 . Thus, We need to pick an inter- 
polation function. In literature several practical forms 
have been used, such as, Bekenstein form fl(x) = (— 1 + 
VI ±4x)/(l + VI + 4z) 0, simple for m p,(x) = x/{l + 

U. Re- 



_4x)/(i + VT 

|29j . and standard form jl(x) 



x) [zy|J, and standard form [x{x) = x/yl + x 
ccntly, a theoretical form from quantum effect in acceler- 
ation universe is proposed fi(x) = (V4x 2 + 1 — l)/2x (30| . 
All these forms can be put into a canonical form [25| • The 
corresponding inverted canonical interpolation function 
is 
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TABLE I. Data of the four selected quasar lens with time de- 
lay measurement. The unit of 0± is arsecond. H band effective 
radius [H is used for HE 2149-2745 and FBQ J0951+2635, and 
V band for SDSS J0946+1835 [|. Ellipticity e = 1-b/a are all 
measured in R band |J l32l , ISH . 



Name 


Zl 


Zs 


e'l 


e'L 


n e 


Ai(days) 


HE 2149-2745 


0.495 


2.03 


1.35 


0.34 0.50 0.5 


103±12 


FBQ J0951+2635 


0.240 


1.25 


0.88 


0.22 


0.17 0.25 


16±2 


SBS 0909+532 


0.830 


1.38 


0.76 


0.42 


1.58 




SDSS J0946+1835 


0.388 


4.81 


2.43 


0.64 


2.59 0.16 


142.7±6.8 



In the following calculations we take the Bekenstein form, 
i.e., (a, 77) = (1,0). 



Although hundreds of quasar lensing have been found, 
only a few have been measured with gravitational time 
delay. This method has been proven difficult because 
the amplitude of quasar variability is quite small, and 
clear and simple modeled relative image is hard to iden- 
tify Up to now there are 19 strong lenses cases with 
measured time delay [3, @l- I n this paper, we consider 
spherical lens only. From the nineteen candidates, we 
select elliptical galaxy lensing systems with double im- 
ages. Only 6 systems satisfy our criteria, they arc HE 
2149-275, FBQ J0951, SBS 0909+532, SDSS0946+1835 
, SDSS J1650+4251,andHE 1104-1805. The rest includes 
clusters, spiral galaxies, multiple images or multiple lens, 
which require analysis beyond the scope of this paper and 
involve more uncertainty. Moreover, the lens galaxy in 
SDSS J1650+4251 is very dark. It does not have a reli- 
able effective radius. And the uncertainty in time delay 
in HE 1104-1805 is too large. That brings our sample to 
a total of four systems. Table Q] lists some properties of 
the four selected systems. 

The luminosity profile of the lens elliptical galaxies sat- 
isfies the de Vaucouleurs' profile. We, therefore, adopt 
the Hernquist mass density profile [34|. The Newto- 
nian gravitational potential and acceleration of Hernquist 
model are $ N = -GM/(g+g h ) and d$ N /dg = GM/(g+ 
Ph) 2 - The Hernquist radius g^ is 0.551 times the effective 
radius (half-light radius). The MONDian acceleration is 
d$/dg = 0(x^)d^/dg (where in = \d^/dg\/ao). 

Different cosmological models give different angular 
distance for the same rcdshift. As MOND is aiming at 
doing away with dark matter, thus, a cosmological model 
consistent with MOND is needed. MOND has been tradi- 
tionally criticized that it can not form large-scale struc- 
ture. Basically, the criticism originated from GR with 
baryons only. However, TeVe S contain additional fields 
which can form large-scale structure compatible with ob- 
servations [35J. As for CMB, TeVeS needs 2eV massive 
neutrinos [36[ which is treated as non-relativistic par- 
ticle, (Qb^Oa) = (0.05,0.17,0.78). We deem that 
Skordis et al. [36| gave a cosmological model consistent 
with MOND. The discrepancy between the two models 
is small in our sample. The difference in D^s / '{D^Ds) 
is within 0.4%. 
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TABLE II. The evaluated mass of the lens and Hubble constant. 
The unit of mass is 10 10 solar mass and the unit of Hubble constant 
is km s — 1 Mpc -1 . x = a/ao is a measure of the MONDian regime 
at the closest approach po> and x± correspond to 8±. 





Mass 




H 








Name 


GR 


MOND 


GR 


MOND 


X- 


z+ 


HE 2149-2745 


no r\25.8 


i6.iil:l 


*7<-) r>81.7 
f Z.Z64.7 


57.6§?:8 


11.3 


2.3 


FBQ J0951+263 


2.9|:1 


2.3|;f 


93.2|^g 5 


79.369i| 


28.5 


4.0 


SBS 0909+532 


-7-7 o78.9 
' ' -^58.3 


cr q56.8 

00.044.3 


84.5g2 1 7 9 


70 9 95 
'U.^68.5 


8.7 


5.9 


SDSS J0946+1835 


95.3go;7 


64.l|f;| 


84.9|fi 


67.og;S 


7.8 


2.8 



In this treatment, we don't have any free parameter. 
Only mass and Hubble constant are unknown variable 
to be solved. Flat rotation curve of spiral galaxies and 
Tully-Fisher relation in gas rich galaxies give a consistent 
value of the acceleration constant cto = 1.21 x 10~ 10 m 
s -2 [H, 14, 2Sj . As we select lensing galaxies that satisfy 
Hcrnquist's model and the Hcrnquist radius would be 
determined by photometric measurement of the effective 
radius, the only unknown in the mass model is the total 
galactic mass. On the part of gravitational lensing, the 
angular position and the redshift (of the lens and the 
two images) and the time delay difference between the 
two images are measured. The only unknown left in the 
lensing part is the Hubble constant. We can evaluate two 
unknowns by solving the lens equation and time delay 
equation simultaneously. The result is summarized in 
Table |Hj The smaller numbers in the first four columns 
are results taken into account of corresponding upper and 
lower uncertainties. 

x = a/ao hi the last two columns of Table [TTI is the ratio 
of the acceleration at the closest approach to the MOND 
acceleration constant. Recall that x 3> 1 is the Newto- 
nian regime and x -C 1 is the deep MOND regime. On the 
one hand, x < 1, MOND gives less mass than GR. The 
excess mass in GR can be interpreted as missing mass or 
dark matter. On the other hand, if x ^s> 1, MOND and 
GR should give similar mass. As shown in Table HH most 
x is of the order of unity (intermediate MOND regime). 
As expected the mass and Hubble constant from MOND 
are smaller than the value computed in GR without dark 
matter. Once again, as in other galactic scale dynami- 
cal phenomena, MOND gives a consistent picture of ex- 
plaining the observed excess acceleration in gravitational 
lensing including time delay phenomenon. 

The other source of uncertainty comes from the choice 
of interpolation function. In Newtonian regime x 1 
or the deep MOND regime x -C 1, different interpola- 
tion functions should give the same result. However, our 
sample lies in the intermediate MOND regime. Table [TTI 
shows the result from Bekenstein form. Other interpola- 
tion functions would give somewhat different result. In 
fact, major uncertainty comes from observation, in par- 
ticular, the time delay measurements. 



In summary, this work is a first attempt to use MOND 
to interpret data from gravitational time delay The Hub- 
ble constant evaluated from this study is consistent with 
literature (see Table |Hj . A simple average of these four 
cases give H n = 68.5 km s _1 Mpc _1 . When compare 
with GR (without dark matter), the evaluated mass of 
the lens is 28% to 49% smaller than those from GR, and 
the Hubble constant is 18% to 27% smaller than GR. We 
stress that the ellipticity of HE 2149-2745 is too large to 
consider it as a spherical potential. We should interpret 
its result with caution. Among the four systems we have 
studied, the data from SDSS0946+1835 is the best as 
it has the least uncertainty in time delay, which is about 
5%. The uncertainty it contributes to mass is about 3.7% 
and Hubble constant about 5%. Perhaps the longer the 
time delay, the smaller the measurement uncertainty will 
be. We look forward to more data from lens with long 
time delay. Gravitational lensing promises to provide a 
testing ground for modified gravity. 
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